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(Dated: December 7,2006) 

We report on a search for second generation leptoquarks (LQ 2 ) which decay into a muon plus quark 
in pp collisions at a center-of-mass energy of y/s — 1.96 TeV in the DO detector using an integrated 
luminosity of about 300 pb -1 . No evidence for a leptoquark signal is observed and an upper bound 
on the product of the cross section for single leptoquark production times branching fraction /3 
into a quark and a muon was determined for second generation scalar leptoquarks as a function 
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of the leptoquark mass. This result has been combined with a previously published DO search 
for leptoquark pair production to obtain leptoquark mass limits as a function of the leptoquark- 
muon-quark coupling, A. Assuming A = 1, lower limits on the mass of a second generation scalar 
leptoquark coupling to a it quark and a muon are ?tilq 2 > 274 GeV and 771lq 2 > 226 GeV for [3 = 1 
and j3 — 1/2, respectively. 



The observed symmetry in the spectrum of elementary 
particles between leptons and quarks motivates the exis- 
tence of leptoquarks [l[ . Leptoquarks are bosons carrying 
both quark and lepton quantum numbers and fractional 
electric charge. Leptoquarks could in principle decay into 
any combination of a lepton and a quark that carry the 
correct charge. Experimental limits on lepton number vi- 
olation, on flavor-changing neutral currents, and on pro- 
ton decay, however, lead to the assumption that there 
would be three different generations of leptoquarks. Each 
of these leptoquark generations couples to only one quark 
and one lepton family and, therefore, individually con- 
serves the family lepton numbers 0]. In this Letter, sec- 
ond generation leptoquarks refer to leptoquarks coupling 
to muons. Since there is no explicit connection between 
a given lepton generation with any of the three quark 
generations in the standard model, the second genera- 
tion leptoquark that couples to muons could couple to a 
quark from any one of the three generations. 

Figure [1] shows mechanisms for leptoquark production 
and decay in pp collisions. Leptoquarks can be either pair 
produced via the strong interaction or single leptoquark 
can be produced in association with a lepton. The cross 
section for single leptoquark production depends on the 
a priori unknown leptoquark-lepton-quark coupling A. In 
pp collisions, the production cross section for single lep- 
toquarks coupling to up and down quarks is significantly 
larger than that for single leptoquarks coupling to sec- 
ond generation quarks, and for the search described in 
this Letter, we only considered this scenario. For other 
quark flavors, the inclusion of single leptoquark produc- 
tion would not improve the sensitivity from the pair pro- 
duction search even for large couplings. 

This search was performed assuming both leptoquark 
pair and single production contribute to the expected 
signal. Therefore both the final state with two jets and 
two muons and the final state with two muons and one 
jet were considered. The former has been studied in 
previously published analyses of leptoquark pair produc- 
tion H H| ■ In addition, limits are given if one assumes 
that only single leptoquark production contributes to the 
expected signal. The cross section limit for this scenario 
can be interpreted as limit on a final state containing two 
energetic muons and a high Et jet. The inclusion of sin- 
gle leptoquark production in a Tevatron search has been 
previously discussed in Ref. [f|. 

The DO detector [6j consists of several layered ele- 
ments. First is a magnetic central tracking system which 
is comprised of a silicon microstrip tracker (SMT) and 
a central fiber tracker (CFT), both located within a 
2 T superconducting solenoidal magnet. Jets are re- 
constructed from energy depositions in the three liquid- 
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FIG. 1: Examples of leading-order Feynman graphs for pair- 
production (left) and single-leptoquark production (right) of 
leptoquarks. 



argon/uranium calorimeters: a central section (CC) cov- 
ering pseudorapidities, rj = — ln[tan(0/2)], where 9 is the 
polar angle with respect to the proton beam direction, up 
to w 1, and two endcap calorimeters (EC) extending 
coverage to \rj\ w 4, all housed in separate cryostats 0]. 
Scintillators between the CC and EC cryostats provide 
sampling of developing showers at 1.1 < \f]\ < 1.4. A 
muon system Q resides beyond the calorimetry and con- 
sists of a layer of proportional wire tracking detectors 
and scintillation trigger counters before 1.8 T toroids, fol- 
lowed by two similar layers after the toroids. The muon 
system is used for triggering and identifying muons. The 
muon momenta are measured from the curvature of the 
muon tracks in the central tracking system. 

The data used in this analysis were collected between 
August 2002 and July 2004, corresponding to an inte- 
grated luminosity of 294 ± 19 pb -1 . Only events which 
pass single- or di-muon triggers were considered. At the 
first trigger level, a muon was triggered by a coincidence 
of hits in at least two of the three scintillator layers of 
the muon system within a time window consistent with 
muons coming from the interaction point. At the second 
trigger level, a muon track was identified from the hits in 
the proportional wire tracking detectors and the scintilla- 
tors of the muon system. The overall trigger efficiency for 
fij + and /i + (ij events fulfilling the selection criteria 
of this analysis was measured to be (89 ± 3)%. 

Muons in the region |?7 < 1.9 were reconstructed from 
hits in the three layers of the muon system which could 
be matched to isolated tracks in the central tracking sys- 
tem. Cosmic muon events were rejected by cuts on the 
timing in the muon scintillators and by removing back- 
to-back muons. Jets were reconstructed using the itera- 
tive midpoint cone algorithm Q with a cone size of 0.5. 
The jet energies were calibrated as a function of the jet 
transverse energy (Et) and r\ by imposing transverse en- 
ergy balance in photon-plus-jet events. Only jets which 
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were well-contained within the detector were considered 
by requiring \rj\ < 2.4. 

For this search, the background is dominated by Drell- 
Yan production and Z decays: Z/j* — * fifi (+jets) 
(Z/DY). Additional backgrounds coming from QCD mul- 
tijet production and from W+jets events (with at least 
one reconstructed muon originating not from the hard 
scattering) were estimated and found to be negligi- 
ble. To evaluate the contribution from the Z /DY back- 
ground, samples of Monte Carlo (MC) events were gen- 
erated with pythia (Version 6.202) [l0| . Samples of tt 
(m t = 175 GeV) and WW production were also gener- 
ated with pythia. The signal efficiencies were calcu- 
lated using pythia samples of LQ 2 + /x — » fj,j + fj, and 
LQ 2 LQ, — > pij + fij MC events for leptoquark masses 
(m L Q 2 ) from 140 to 280 GeV in steps of 20 GeV. All MC 
events were processed u sing a full simulation of the DO de- 
tector based on GEANT [1JJ and the standard event recon- 
struction. Differences in the trigger and reconstruction 
efficiencies between data and Monte Carlo were taken 
into account using proper weightings of the MC events. 

The search for leptoquarks required two muons with 
transverse momenta pt > 15 GeV and either one or two 
jets with E T > 25 GeV. To reduce the Z/DY background 
at high dimuon mass due to occasionally poorly recon- 
structed muon tracks, advantage was taken of the fact 
that no or little missing transverse energy is expected 
in either signal or Z /DY events. The missing transverse 
energy was determined from the transverse energy imbal- 
ance of all muons and jets (Et > 20 GeV) in the event, 
and the momentum of the muon opposite to the direc- 
tion of the missing transverse energy (i.e. with the larger 
azimuthal angle relative to the direction of the missing 
transverse energy) was corrected such that the missing 
transverse energy parallel to the muon vanished. This 
degraded the muon momentum resolution and shifted the 
dimuon mass to lower values in both data and MC. How- 
ever, this correction allowed the suppression of poorly 
reconstructed Z boson events shifted into the high mass 
region where the search for leptoquarks was taking place. 

To create statistically independent signal bins, events 
were first classified as either leptoquark pair or single 
leptoquark candidates. Events were classified as lepto- 
quark pair candidates if they contained two jets with 
Et > 25 GeV, had a dimuon mass M PM > 105 GeV (to 
remove Z boson events), and fulfilled the requirement 
S = S T / GeV - 0.003 x (M w / GeV - 250) 2 > 250, with 
St defined as the sum the sum of the absolute values 
of the transverse energies of the two jets and the trans- 
verse momenta of the two muons forming the fij + fij 
system. Events not classified as leptoquark pair candi- 
dates were classified as single leptoquark candidates if 
they contained at least one jet with Et > 50 GeV, had a 
di-muon mass M m > 110 GeV and fulfilled the require- 
ment E = (M^/GeV - 110) x (E T /GcV - 50) > 500 
(see Fig. [5]). The optimum choice of variables and cut 
values has been determined to optimize the sensitivity 
to the signal. These selections have been cross-checked 




[GeV] 




40 60 80 100 120 140 160 180 200 



E J T [GeV] 

FIG. 2: Distributions of M MM vs. the highest jet E 3 T for back- 
ground (upper plot) and signal and data (lower plot). The 
curved lines show the edges of the signal bins (see text for 
definition). The bottom plot shows as small dots the ex- 
pected signal distribution from single leptoquark production 
(tilq 2 = 260 GeV) . The bottom plot also shows the observed 
data for events not classified as leptoquark pair candidates 
(small black boxes) and classified as leptoquark pair candi- 
dates (open boxes). 



with a neural net optimization, which gave similar re- 
sults. Eleven events were either classified as leptoquark 
pair candidates or single leptoquark candidates while 
6.6 ± 0.5 (stat) ± 1.1 (syst) are expected from standard 
model background. A small excess of data over back- 
ground was observed. The probability that 6.6 ±1.2 ex- 
pected events fluctuate up to 11 observed events is 9.2%. 

Candidate events were arranged in bins of increasing 
signal to background ratio. For leptoquark pair candi- 
dates, bin boundaries of S = 320 and 390 are used to 
define bins PI, P2 and P3 For single leptoquark can- 
didates, boundaries of E = 1000, 2500, 5000 were used to 
define signal bins SI, S2, S3 and S4 (Fig.©. 

Table Q] summarizes the efficiency of the single lepto- 
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TABLE I: Signal efficiency (e a ingie) for selecting single leptoquarks for (3 = 1, number of data events (JVdata), and number 
of predicted background events (N^^ d ). The first uncertainty on N^^ d is due to limited Monte Carlo statistics, the second 
denotes the systematic uncertainty. The first two lines indicate the total number of events after the initial event selection while 
the other lines indicate the numbers for the individual bins of the leptoquark pair candidates (P1-P3) and single leptoquark 
candidates (S1-S4) as described in the text. 



Cut 


£ single 

(m L q 2 = 200 GeV) 


£ single 

(m LQ2 = 240 GeV) 


■Wdata 




Mfifj, > UOGeV 


0.145 ±0.013 


0.176 ±0.015 


43 


44.75 ± 1.74 ±6.13 


F, j > 95 dpV 

J-J rp / £iKJ VJC V 










Mfi/j, > UOGeV 


0.122 ±0.012 


0.158 ±0.014 


20 


13.41 ± 0.92 ± 1.57 


E J T > 50 GeV 










PI 


0.011 ±0.002 


0.015 ±0.002 


2 


0.96 ±0.25 ±0.28 


P2 


0.006 ±0.001 


0.011 ±0.002 


2 


0.39 ±0.10 ±0.11 


P3 


0.006 ±0.001 


0.012 ±0.002 





0.27 ±0.10 ±0.08 


SI 


0.018 ±0.002 


0.014 ±0.002 


2 


2.01 ±0.33 ±0.57 


S2 


0.028 ± 0.003 


0.030 ± 0.003 


1 


1.61 ±0.27 ±0.44 


S3 


0.016 ±0.002 


0.029 ± 0.003 


3 


0.87 ±0.17 ±0.29 


S4 


0.015 ±0.002 


0.029 ± 0.004 


1 


0.44 ± 0.08 ± 0.06 


Signal bins (Pl-P3,Sf-S4) 


0.100 ±0.010 


0.140 ±0.013 


11 


6.55 ±0.53 ±1.08 



quark selection for two leptoquark masses as well as the 
number of expected background events and the distribu- 
tion of the data in the three pair and four single lepto- 
quark bins. 

The dominant uncertainties on the predicted number 
of background events are Monte Carlo statistics, vary- 
ing between 7% and 25% for the seven signal bins, the 
jet-energy scale uncertainty [(2 - 12)%], and the uncer- 
tainty on the shape of the jet transverse energy distribu- 
tion of Z/DY events [(20 - 30)%]. The latter has been 
estimated by a comparison of the PYTHIA [l(| simula- 
tion with Monte Carlo events generated with the ALP- 
GEN [l2j event generator. For leptoquark pair candidates, 
the uncertainty due to the two jet requirement was esti- 
mated to be 16% Q. In addition, the following sources 
of systematic uncertainties were taken into account: lu- 
minosity (6.5%), theoretical cross section of the Z/DY 
processes (3.6%), and muon triggering and identification 
(5%). The systematic uncertainties, added in quadra- 
ture, are shown in Table [IJ 

The systematic uncertainties on the signal efficiencies 
arise from the jet-energy scale uncertainty [(3 - 17)%], 
muon triggering and identification (5%), limited Monte 
Carlo statistics [(4 - 14)%, uncorrelated between bins], 
and parton distribution function uncertainty (2%). 

Limits on the cross section for single leptoquark pro- 
duction were calculated from the observed and expected 
events in the seven signal bins (P1-P3 and S1-S4). Three 
different scenarios were considered: (a) the only contri- 
bution to the signal region came from standard model 
background and from single leptoquark production, (b) 
contributions to the signal region came from standard 
model background and single leptoquark production plus 
leptoquark pair production corresponding to the nominal 
leptoquark pair cross section [13| with (3 = BR(LQ — > 
fxq) = 1/2, and (c) as (b) but with (3=1. For scenar- 
ios (b) and (c) leptoquark pair events in the signal bins 
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280 



M LQ [GeV] 



FIG. 3: Cross section upper limits at 95% C.L. for the pro- 
duction of single leptoquarks for the three scenarios (a) no 
contribution from leptoquark pair production, (b) pairs con- 
tribute with a (tx/3 2 corresponding to (3 — 1/2, and (c) pairs 
contribute with a ax/3 2 corresponding to (3 = 1. The solid line 
is the observed limit and the dashed line the average expected 
limit assuming that no signal is present. Also indicated are 
the predicted single leptoquark production cross sections for 
A = 1, p = 1 and A = 1, j3 = 1/2. The shaded region is 
the variation of the cross section using renormalization and 
factorization scales of 2 x ?tilq 2 and 0.5 x itilq 2 , respectively. 



P1-P3 are treated exactly the same way as in 3]. The 
analysis described above can therefore be considered as a 
combination of the search for singly produced leptoquark 
with the leptoquark pair analysis published in [3j. 

These calculations were performed assuming a flat 
prior and Gaussian errors as described in Ref. LJ] with 
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— / / 
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u 180 200 220 240 260 280 

M LQ [GeV] 

FIG. 4: Upper limits on A 2 for the three scenarios: (a) no 
contribution from leptoquark pairs events and /3 = 1, (b) 
leptoquark pairs contribute with a a x 1 corresponding to 
j3 = 1/2, and (c) leptoquark pairs contribute with a a x /3 2 
corresponding to P = 1. The solid line is the observed and 
the dashed line is the expected limit. The regions above the 
solid lines are excluded at 95% C.L. 

the correlations between the systematic errors taken into 
account. The observed limit was calculated using the 
confidence level ratio [3] CLs = CLs+b/CLb, where 
CLs+b is the confidence level for the signal plus back- 
ground hypothesis, and CLb is the confidence level for 
the background only. 

Figure [3] shows the 95% C.L. exclusion limits for the 
production cross section times branching fraction for sin- 
gle leptoquarks as functions of the leptoquark mass. The 
predicted cross section for single leptoquarks depends 
both on the coupling strength A and on the leptoquark 
mass. Each point in Fig.[3]therefore corresponds to a spe- 
cific value of A and branching fraction. In scenarios (b) 
and (c), leptoquark masses are excluded independently 
of the single leptoquark cross section if they are already 
excluded because of the leptoquark pair production re- 
sults. The average expected limits are better than the 
observed ones because of the small excess of data over 
background. 

Also indicated are the predicted leading order cross 



sections for single leptoquarks [l5[ for A = 1 with (3 = 1 
and (3 = 1/2. For the production cross sections, it was 
assumed that the leptoquark couples to the u quark and 
a muon. From the intersection of the cross section limit 
with the lower limit of the predicted cross section, which 
were derived by using renormalization and factorization 
scales of 2 x tolq 2 , respectively, exclusion limits on the 
leptoquark mass as a function of the coupling strength 
A can be estimated. Figure [4] shows the 95% C.L. ex- 
clusion regions as functions of the leptoquark mass and 
A 2 . For scenario (b), the reduction of expected events 
due to (3 = 1/2 applies to both single and pair produc- 
tion of leptoquarks. The intersection of the curves with 
A 2 = yields the result when only leptoquark pairs are 
considered. 

The production of a leptoquark-like particle in associ- 
ation with a muon which decays into a jet and a muon 
can be excluded for cross section times branching frac- 
tions a x (3 > 0.40 pb for a particle mass of 200 GeV 
and a x [3 > 0.23 pb for a particle mass of 260 GeV. 
The DO Run II result for scalar leptoquarks and (3 = 1 
°f m LQ 2 (A 2 < 1) > 247 GeV, which only considered 
leptoquark pair production 3], is improved for an as- 
sumned leptoquark coupling to a u quark and a muon 
of A 2 = 1 to m LQ . 2 (A 2 = 1) > 274 GeV. For f3 = 1/2, 
the improvement is from mi J Q 2 (X 2 « 1) > 190 GeV to 
™lq 2 (A 2 = 1) > 226 GeV. For A 2 = 0.1 the observed 
limits show no improvement while the expect limits in- 
crease by about 7 GeV. This analysis is the first search 
for single leptoquark production in pp collisions. 
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